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Abstract 

We investigate the nonlinear dynamics underlying the evolution of a 2-D nanoscale 
ferromagnetic film with uniaxial anisotropy in the presence of perpendicular pump- 
ing. Considering the associated Landau-Lifshitz spin evolution equation with Gilbert 
damping together with Maxwell equation for the demagnetization field, we study 
the dynamics in terms of the stereographic variable. We identify several new fixed 
points for suitable choice of external field in a rotating frame of reference. In par- 
ticular, we identify explicit equatorial and related fixed points of the spin vector in 
the plane transverse to the anisotropy axis when the pumping frequency coincides 
with the amplitude of the static parallel field. We then study the linear stability 
of these novel fixed points under homogeneous and spin wave perturbations and 
obtain a generalized Suhl's instability criterion, giving the condition for exponential 
growth of P-modes under spin wave perturbations. Two parameter phase diagrams 
(in terms of amplitudes of static parallel and oscillatory perpendicular magnetic 
fields) for stability are obtained, which differ qualitatively from those for the con- 
ventional ferromagnetic resonance near thermal equilibrium and are amenable to 
experimental tests. 
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1 Introduction 



Interest in ferromagnetic resonance (FMR) has soared in recent times due to 
advances in fabricating nanostructures. This imphes prospects for several new 
experiments to study possible absorption phenomena, and patterns that may 
form owing to instabilities, in ferromagnetic films. While other magnetic reso- 
nance counterparts such as nuclear magnetic resonance (NMR), electron para- 
magnetic resonance(EPR), electron spin resonance, etc., have found immense 
technological applications, including crystal structure determination and med- 
ical diagnostics, FMR has remained a more complex phenomenon. Some of the 
main features of FMR are a) large magnetization, and hence large susceptibil- 
ity, b) large demagnetization field due to strong magnetization, which is also 
influenced by the shape of the specimen, and c) resonance excitations breaking 
into spin wave modes, that make spin reversal more complex. 

Spin-wave instabilities were first observed independently by Damon [1] and by 
Bloembergen and Wang [2] as noisy anomalous absorption which abruptly sets 
in at a certain threshold power as the resonant microwave field is increased. 
Suhl remarked that this phenomenon "bears a certain resemblance to the 
turbulent state in fluid mechanics" [3] . In fact the instabilities in this case 
were first explored in [3] and are referred to as Suhl instability. These were 
subsequently verified experimentally [4]. 

Under a growing attention to deterministic chaos and nonlinear dynamics, 
there occurred a renaissance on studies involving high-power magnetic reso- 
nance in the 1980s. High-resolution experiments were carried out for spin-wave 
nonlinear dynamics in the yttrium iron garnet film and sphere in parallel 
and perpendicular pumpings [5,6,7,8]. For certain high powers beyond the 
Suhl threshold, interaction between excited spin-wave modes lead to various 
dynamical phenomena including auto-oscillations, period- doubling cascades, 
quasi-pcriodicity, and chaos. Also observed were irregular relaxation oscilla- 
tions and abrupt transitions to wide-band turbulence, beyond the Suhl thresh- 
old [5,6,7,8], which were explained by using Zakharov et al's S-theory [9]. Since 
then the studies on spin wave instabilities have acquired renewed interests up 
to now [10,11,12]. 

Despite these pioneering works, including that of Suhl[3], the investigations 
have been limited to the instability around fixed points that correspond to 
magnetization parallel to the anisotropy axis. This is due to the following fact: 
In contrast to the case of NMR and EPR, macroscopic ferromagnets in FMR 
have large frozen magnetization. Under a static magnetic field, therefore, any 
available pumping field cannot freely rotate such strong magnetization at the 
cost of the large stabilizing energy (Zeeman energy). In nanoscale ferromag- 
nets, however, it is possible to rotate the saturation magnetization far from 
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the anisotropy axis. The magnetization can even be driven to the equatorial 
plane perpendicular to the anisotropy axis. Therefore it is timely to analyze 
the nonlinear dynamics of spin waves in these new non-equilibrium states. 

It may also be noted that a number of investigations exist in the literature 
on the dynamics of higher dimensional Heisenberg ferromagnetic spin systems 
[13,14,15] corresponding to isotropic (pure exchange interaction), anisotropic, 
external field and other interactions [13,14,15,16,17]. However, to the knowl- 
edge of the authors, concerning the systems coupled with Maxwell equation 
for the demagnetization field as considered as in this paper, there exist very 
few studies. 

In this article, we investigate a 2-D ferromagnetic system with uniaxial anisotropy 
in a background of alternating magnetic field transverse to, and a static mag- 
netic field parallel to, the anisotropy axis. The axis of anisotropy is chosen 
to be arbitrary. We also include the demagnetization field due to the spatial 
inhomogeneity of the magnetization vector, which is seen to play a crucial 
role. Fixed points - P-modes - of such a system have been identified earlier, 
and their stability against both homogeneous perturbations and spin-wave 
excitations analyzed [18,19]. However, we specifically analyze equatorial and 
other fixed points, that have not been identified earlier, wherein the average 
magnetization vector lies in the plane transverse to the anisotropy axis and 
exhibit a more complex dynamics. We further obtain a criterion for instability 
under spin wave excitations, thus generalizing the Suhl's instability criterion. 

The strategy we follow in our analysis is as follows. We project the unit spin 
vector m(r,t) stereographically onto a complex plane ip{r,t) and deduce the 
equation of motion in terms of the stereographic variable. Now going over 
to a natural rotating frame of reference, one is able to identify the defining 
equations for the fixed points or the so called P-modes. Of all the possible equi- 
librium points, the equatorial and related fixed points are of special interest as 
their expressions can be explicitly obtained. In fact we find that there exists 
four such equilibrium points. Then we investigate their linear stability nature 
(i) under spatially homogeneous perturbations and (ii) under more general 
spin wave perturbation in order to identify their local and global stability, 
and obtain conditions for Suhl's instability as a function of experimentally 
measurable parameters such as the amplitudes of oscillatory perpendicular 
and static parallel external magnetic fields. The results give clear criteria for 
experimental realization of the predicted results. 

The plan of the paper is as follows. In Section 2, we introduce the model spin 
Hamiltonian for the ferromagnetic film in the presence of the external field and 
an intrinsic demagnetization field, and write down the Landau-Lifshitz (LL) 
equation for the spin field with the Gilbert damping term included. We further 
introduce the stereographic variable, and rewrite the LL equation in terms of 
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the new variable. In Section 3, we identify fixed points of the LL equation, 
specifically the equatorial and other related fixed points, and analyze their 
linear stability under homogeneous perturbations. In Section 4, we study the 
stability of these fixed points under spin wave excitations in terms of a period 
map and generalize Suhl's instability criterion. We conclude with a summary 
of results in Section 5. 



2 Model Hamiltonian and LL equation 



We consider a 2-D ferromagnetic film with uniaxial anisotropy and an ap- 
plied oscillating magnetic field in the direction transverse to the anisotropy 
axis. Such a system can be described by a Hamiltonian for the unit spin field 
m(r, t) = {mi(r, t), m2(r, t), m3(r, t)}; = 1 as 

H — Hgxchange ~l~ H^ppHgf^ -\- Ha^iigofj-gpy -\- Hd^^agYietizatiom (1) 



where 

Hexchange = / D(Vnif dh, (2) 

Happlied = - j Ba • ih (pT, (3) 

Hanisotropy — J KUl^^ d T, (4) 

-J • m d^r, (5) 



H demagnetization 



and the demagnetizing field satisfies the Maxwell's equation, 

V • = -47rV • m; V x = 0. (6) 



In the above, we have assumed the average magnetization to be of constant 
magnitude in space and time, and hence conveniently expressed the Hamilto- 
nian (1) in terms of the normalized unit average magnetization ih. Equation 
(6) gives the constitutive relation connecting the demagnetization field with 
the average local magnetization, assuming absence of any free current in the 
film. In equations (2-5), D is the exchange constant, is the applied magnetic 
field, K is the anisotropy parameter and my is the component of magnetization 
in the direction of anisotropy. 
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Time evolution of the spin field is governed by the Hamilton's equation 

r IT 

— = {m,i/}. (7) 



Using the spin Poisson bracket relations 

{mi{r, t),mj{r', t')} = eijkmk5{Y -r',t- t'), (8) 

it is straight forward to show that the dynamics is governed by the LL equation[20], 
m = -m X Feff, (9) 

where Fg// is the effective field given by 

Fe// = i:)V^m + B„ + «;m||e|| +H^. (10) 

In equation (10), is the two dimensional Laplacian = d'^ / dx'^ + d"^ / dy'^ . 
The direction of anisotropy ey has been chosen to be arbitrary in the internal 
space of the spin field. Further, writing ey = {0, 0, 1}, we choose the applied 
magnetic field to be of the form 

Ba = ha±{cosu!t, sina;^, 0} + ha\\e\\ = ha±{cosu!tei + smu!te2) + ^a||e||, (11) 



where Ci and 62 are two unit orthonormal vectors in the plane transverse to the 
anisotropy axis. Notice that the orientation of the anisotropy axis is specified 
as in equation (11) in the spin space, but is still arbitrary in the physical 

space, where the plane of the ferromagnetic film is assumed to be in the x — y 
direction. Without loss of generality one can choose my = m^. Introducing a 
phenomenological Gilbert damping term (proportional to a), equation (9) can 
be modified to the form 

m — am x m = — m x Fg//. (12) 



It is this equation in 2-D coupled with the Maxwell equation which we shall 
investigate in this paper. 

It proves convenient for later analysis to express the equation in terms of the 
complex stereographic variable (see Figure 1 for details) 

_ mi +im2 , , 

1 + ms 
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or equivalently by the transformation 



(14) 



Correspondingly, the terms in the Hamiltonian in equation (1) become 

H exchange ~ 



applied 



-I 



h. 



1 + 



'a±- 



H anisotropy 



K- 



l + IVP 



+ ha\\ 



(i+iV'P) 



d'r, 



H, 



demagnetization 



1+ |^|2)2 



1 + 



+ K 



(i-IV^P) 



(i + IV'P) 



(15) 
(16) 
(17) 
(18) 



where we have used equation (11), and represented the demagnetization field 

as Urn = {^1, ^2, hm\\} with Hjn = /ll + i/l2- 

On account of the spin Poisson bracket relations in equation (8), ^ obeys the 
Poisson bracket relation, 

{^|;iv,t),r{r',t')}^^Sir-r',t-t'). (19) 



For the complex scalar field ip, one can write down the evolution equation from 
the Hamiltonian (1), equivalent to equation (9). Then, including the damping 
term, the evolution equation for ip equivalent to equation (12) reads as [21] 



i{l - ia)ip = D{V^ij; - 2 



1 + I^P 



K- 



(20) 



-hm\l^ + ^{Hr^-iP^H*J. 



The form of the applied transverse field in equation (20) suggests that it is 
convenient to make a transformation 

ip = ip{r, t) exp[iujt], (21) 
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and rewrite equation (20) in terms of -0- This amounts to moving to a rotating 
coordinate frame in the spin space, rotating about the anisotropy axis ey with 
angular frequency uo. If m were expressed in spherical polar coordinates in the 
rotating frame as m = {sin 9 cos 0, sin 9 sin 0, cos 6'}, then ip = tan(6'/2) exp 
With this change, the equation (20) becomes 



i(l - ia)-^ = L>(V^V' - r ^ y ) - {hai\ -u; + iau; + ^)^ (22) 

It must be remembered that hm\\, Hm and in equation (22) are intrinsically 
related to -0 through the relation (6). Due to the second equation in (6), we 
can write 

= -V(^, (23) 



by introducing an auxiliary potential 99(r, t). Thus, can be found by solving 
the Poisson equation V^</? = — 47rV • ih. Consequently, the boundary condi- 
tion and shape of the boundary influence the actual form of </?, and hence 
Hm- In the following we shall consider only spatially homogeneous spin fields 
and small amplitude spin waves. In this case, boundary effects can be ignored 
if one assumes the spin wave to be confined to a localized region, and that 
they dissipate close to the boundary. In the case of large amplitude perturba- 
tions, coupling between different modes become nontrivial. This case will be 
separately treated elsewhere. 



3 Fixed points (P-modes) and stability 



Spatially homogeneous fixed points of equation (22) correspond to a uniform 
magnetization field that exhibit a periodic motion of frequency a;, irrespective 
of the fact that the governing equations are highly nonlinear. Such homoge- 
neous steady states in the rotating frame are referred to as P-modes. These 
fixed points ipo-, in the moving frame, are obtained from equation (22) by as- 
suming time independence and spatial homogeneity of ■0, i.e., dipo/dt = = 
Wipo- For such a uniform steady field ipo, the uniform demagnetization field can 
be conveniently written as = — A^j_mj^ — A'"||m|| , where m° = {m°, 0}, 
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I.e., 



H^n. = -N^rnl + iml) = (24) 

1 - h/; P 

/iHI = -iV||m|| = -iV||-— (25) 



Using these conditions in equation (22), one obtains after some rearrangement, 
the equation for the fixed points as 



(6„ll + + \ - \k^{1 - V'o) = 0, (26) 



where the quantities ha\\i ba±, fi and n^ff are defined as 

ha\\ - ^ , o '^'^ , Ar AT ^07^ 



Evidently, the fixed points are functions of the apphed field strengths ha_L and 
h\\, and the frequency a;, besides the anisotropy constant k. However, of these 
h\\, hai_ and cu are the tunable parameters. 

From equation (26), it is straightforward to check that the equation satisfied 
by iV^oP is. 

Equation (28) can be rewritten as a polynomial of degree four in u — |'^op as 



- ^li^ + (fe^ii - 26„|| + 1 + + (26^11 _ 2 - 2^2 + (29) 

+(6^11 + + 1 + 2ba\\)u - ^ = 0. 

The quartic equation (29) admits in general zero, two or four real solutions 
depending on the sign and magnitude of the coefficients. These roots can 

be easily obtained for any given set of parameters, and the nature of the 
corresponding fixed points can be analyzed by means of standard theory of 
equilibrium points of two-dimensional dynamical systems. One can establish 
the presence of either stable modes or foci, unstable nodes or foci, or saddles. 
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A detailed analysis of these modes in terms of polar coordinates has been 
made in the past by Bertotti et. al. [18]. 

In this paper, we wish to concentrate on certain physically interesting modes 
not identified earlier, namely those arising out of fixed points along the equator 
or related fixed points, whose form can be explicitly determined. Consequently 
their stability analysis can also be carried out in detail. In particular, on ac- 
count of the anisotropy, the tendency for the spin field is to orient along the 
direction of iey . Naturally, the fixed points along the equator, i.e. , in the 61—62 
plane, are special, being the other extreme, but are available in nonequilibrium 
nano-ferromagnets. Besides, as will be shown below, these fixed points are de- 
generate, and the nature of dynamics around these fixed points is naturally of 
considerable interest. In the next three sections, we identify these fixed points 
and study their stability under homogeneous and spin wave perturbations. In 
Section 4.2, we obtain the generalization of the Suhl's instabihty criterion. 



3. 1 Equatorial and related fixed points 



In this section we identify fixed points along the equator and others related 
to them, and analyze their stability under homogeneous perturbations. 

Along the equator, we have ra°\^{= mg) = 0, or \'ipo\^ = 1- Conversely, |'?/'oP = 1 
implies spins are pointed in the equatorial plane. But from equation (28), fixed 
points along the equator are possible only when ba\\ = 0,or ha\\ = a;. However 
for this choice of the pumping frequency cu or the equivalent amplitude of the 
static parallel field ha\\, other related fixed points also exist. In fact with this 
choice, equation (26) can be solved to obtain four fixed points Pi, i-*2, -P3 and 
P4: 



olm 



tan (j) = —Q] 



oRe 



Pi = P2: = 1; 



P3 : \i^o\' = 



P4 : iV-ol' = 



2(1 + ^2) r 



1 + ^1 



ft2 
1 + 



2(1 + 1]2^ r 



-'a± 



1; 



- 1. 



(30) 
(31) 
(32) 

(33) 



The fixed points are adjusted by tuning the transverse field strength ha±{= 
><-effK±) and the frequency uj{= Kefffl/a) of the applied field, while keeping 
ha\\ — CO. Equations (31)-(33) also impose a limit on the applied transverse 
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field strength, namely, < < 1 + In terms of the original variables 
this condition becomes < hlj_ < '^lff + o?^"^- 



From equation (28) we notice that there are four solutions to j-^oP as seen 
above. Of these four, two fixed points, {P\^P-^ merge (iV'oP = 1) with the 
choice ha\\ = 0, or ha\\ = uj. These fixed points lie along the equator, transverse 
to e\\. The other two related fixed points (P3, P4) in equations (31-33) vary 
between south pole (I'l/'oP = 00) and equator for the range < h^,^ < 1 + fi^. 
The exact orientation will also be decided by the azimuthal angle determined 
by n, as in equation (30). Particularly, for the choice h"^^ = 1 + Q^, i.e., 
/i^j^ = Kg^j + a^uP', the only possible fixed points are along the equator. 



3.2 Spatially homogeneous perturbations 



For a spatially homogeneous perturbation 5'?/'(t) from ijjo, the linear stability of 
ipo is analyzed using equation (22), with the assumption = ipo + tdip, e << 1, 
and restricting to linear terms in 5ip. Correspondingly, the perturbed equation 
reads as 



1 _ U/; |2 . 

i{l - ia)5^P = -neffiin + P^)# + 2n,ffV'il^o{il^M* + V':^) (34) 



1 + l^ol 



-ha±i'oSi', 



where we have defined F = (1 + j'^op)"^. The equation (34) and its complex 
conjugate can be written in a compact form using an appropriate matrix A, 
as 



A T 



(35) 



with 



l—ia 



in 



2T\^o\ 



2iK, 



eff 



1+ia 



*2 



eff 



1+ia 



1-ia ^ 



-ifl + 



2F|^„ 



(36) 



and ^' = {Sip, Sip*)'^. The type of instability of a fixed point is determined by 
the determinant and trace of A. Defining x = (1 — \'4'o\'^)/{^ + l^oP), from 
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equation (34), the determinant and trace are found to be 



Tr(A) = -2a-^(i±^ + ^). (38) 
1 + a"^ 2 a 



In general, the fixed point is a stable node or foci if |A| > and TrA < 0, 
an unstable node or foci in case |A| > and TrA > 0, and a saddle if 
|A| < 0. However, for the equatorial fixed points iV'oP = Ij and consequently 
|A| = 0, since % = 0, thus leading to degeneracy. In order to identify the 
nature of these fixed points, we directly solve equation (22) numerically in the 
neighborhood of the fixed points, assuming spatial inhomogeneity of ip and plot 
the resultant phase portrait in Figure 2. Figure 2(a) shows the phase portrait 
for the choice fe^j^ = for which the only fixed points are along the equator 

(PiandP2), as remarked earlier. Figure 2(b) is for the choice = (1 + Q'^)/2, 
for which there are two additional fixed points at ipo = (2.4, —0.2) (P3) and 
4 = (0.41,-0.03)(P4). 



4 Equatorial and related fixed points and spin-wave instability 



Unlike in EPR and NMR, spin reversal is prevented in ferromagnetic reso- 
nance because of the large frozen magnetization and of the emergence spin 
wave modes coupling to resonance excitations. In this section we generalize 
the study in [3] to the case of equatorial states of ferromagnet with uniaxial 
anisotropy in an arbitrary direction, and a demagnetization field. In Section 
4.1 we concentrate specifically on equatorial fixed points in the presence of 
spin wave excitations. In Section 4.2 we obtain the generahzed Suhl's instabil- 
ity criterion. For this purpose, we consider first a general perturbation 5ip{r, t) 
to the homogeneous steady field ipo, in contrast to the spatially homogeneous 
perturbations of the previous section. The nature of instabilities are identified 
by linearizing equation (22) in 6ip, around ipo- Upon substituting ip = ipo + ^i' 
in equation (22), and confining to first order in 5il^, we find. 



i(l - ia)5ip = DV^(5V' - (haW -uj + iauj + k ^)# (39) 

1 + l^oP 

+2/tr^„(V'o#* + ro^i') 
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Note that in the above equation (39), the dispersion term and the variations in 
the demagnetization field are nontrivial due to the perturbation being spatially 
inhomogeneous. Making now a spatial Fourier decomposition, 

= E = E (40) 



and substituting in equation (39) we obtain the mode equations for with the 
propagation vector q. Due to the Poisson bracket relations in equation (19), 
the mode amplitudes also obey the Poisson bracket relation {aq,a*_^,} — 

The corresponding change in the demagnetization field is obtained using the 
relation 

V • SU^ = -47rV • 6ni. (41) 



Here, 6m is the change in the average local magnetization corresponding to 
Sip. This can be calculated using equation (13), and is given by 

dm = r| {e-'^' - i^y^')d^* + {e'^' - e-*'"*)5^, (42) 

Using equation (40) in equation (42) we can write, 

5m = EZ,(t)e''^■^ (43) 



where 



^q{t) = r|a*_,(e-^'^* - ipy^') + aqie"^' - tPfe''^'), (44) 



^%it). (45) 



Consequently, it can be verified that 

5H^ = -47rE4(q-Z,)e^''-'" (46) 

q Q 
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satisfies the constitutive relations in equation (6). 

From equation (39), using equations (40) and (46), we iiave tiie mode equation 
satisfied by the amphtudes Oq, 



i{l — ia)aq — — Dq^ — (ha\\ — to + iau + k ) 



(47) 



1 + 

where SH^ and Shm\\ are obtained from equation (46). If q(= q/|q|) is ex- 
pressed in the spin frame as q = sin^g(cos0gei + sin 0^62 ) + cos^^ey, then 



sin^„ ag(e'("*-'^') - 4*'e-'("*-^') 
+a* (e~'('^*~'^«^ — '02gi('^t-</'q)^ 

-2 cos 9 q{'tpoa*_g + iplaq 



(48) 



qZ^ 



(49) 



Here 6q and 0g are the polar and azimuthal angles made by the wave vector 
q in the spin frame spanned by Ci, 62 and ey. Substituting equation (48) in 
equation (46) and using the resultant expression for SUm in equation (47), we 
can write 



(Iq + iAqttq + iBqOLq = 0, 



(50) 



^{l-iaY^ - Dq'^ - (h\\ - u + iau + K- (51) 

L ^ 1 + IV'op 

-27rr|2V'oCOS^q - sin^q(e-^('"*-'^«) - ^2giM-0,)-)|2 



Bq = {l- la)-^ 



2KVijl - 27rr( 27/^0 COS 6, 



(52) 



-sin^q(e-^('"*-'^«) - i)y^'^^-'>"''^))' 
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Similarly, 

a*_g — iA*_ga*_g + iB*_qaq = 0. 



(53) 



Now we will analyze these mode equations (50) and (53) in some detail in the 
following subsections. 



4-1 Stability of equatorial and related 



fixed points As we noted in Section 3.1, the fixed points along the equator. 
Pi and P2, as well as the related fixed points P3 and P4 are possible when 
the pumping frequency u coincides with the amplitude of the static parallel 
magnetic field ha\\ ■ The nature of stability of these fixed points can be studied 
by identifying a period map M using the mode equations (50) and (53). In this 
connection it must be remembered that the P-modes have inherent in them 
a time period T = 27r/a; arising from the frequency of the applied transverse 
field. Then the period map M, defined through the equation. 



^ S(T) 



M 



«%(0) 



(54) 



is obtained from equations (50) and (53) after integration through one period 
T. Numerically, the map M in equation (54) is obtained by evolving over a 
period T the two column vectors (1, 0)"^ and (0, 1)"^, taken as the initial vectors 
(ag(0), alg(O))^. The two columns so obtained then form the columns of the 
period map M. Prom equations (30) and (31)-(33), there are effectively two 
tunable parameters, fl and b± or uj and ha±, that also fix the fixed point ipg. 
If 7g-j_ are the eigenvalues of the period map M for a given fixed point ipo, 
the fixed point is unstable for a mode of wave vector q if |7c,±| > 1. Figure 3 
shows unstable regions in the cos 9 — u space as shaded regions for different 
values of the angle 9q: (a) sin^^ = 0, (b) sin^g = 0.6 and (c) sin^^ = 1. Note 
that cos^ here effectively stands for ha± through the relations (31)-(33). In 
Pigure 4, the same regions arc shown in the ha± — ha\\ plane. However, not all 
points on the plane are experimentally accessible if one requires a homogeneous 
background, due to the condition < < n^^f + a^a;^. Note that 6q defines 
the angle between the anisotropy axis and the plane of the ferromagnetic film. 
It should be emphasized that the global stability diagrams in Pigures 3 and 
4 are obtained under the condition that uo = ha\\. We find that the equatorial 
and related fixed points are always stable for the perpendicular field with 
its strength ha± less than the value determined by the pumping frequency, 
but become unstable when ha± is "decreased" below the generahzed Suhl 
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threshold. There exists an interesting singular structure near the vicinity of 
u; = (Figure 3) or of ha\\ = (Figure 4), which is very sensitive to the angle 
6q. These phase diagrams differ qualitatively from those for the conventional 
spin-wave instability near the ferromagnetic ground state. 



4-2 Mode equations and SuhVs instability 



It can be verified from equations (51) and (52) that Ag 



A_q and Bq 



B-q, if 



we note q 



q 



TT- 



■ 9q and 



. + 71. To begin with, we shall first 



suppress the oscillating terms with exp[±i{ujt — 0^)] and exp[±2i(Ljt - 
Then, Aq and Bq in equation (50) are replaced by Aq and Bq given by 



Aq — {1 — ia) ^ — Dq^ — ^/ly — u + iauj + k-^ 



1 - |4| 



-27rr (4 1 4 r cos^ Oq + sin^ ( 1 + | ^ | ^)) 



(55) 



Bq^{l- ia)-^ 2«;rVio - 27rr(4Vio cos^ Oq - 2 sin^ Oqi/jl 



(56) 



Then equation (47) can be written in terms of normal modes, after a Bogoli- 
ubov type linear transformation aq = Xqbq — iXqh*_^^ to a new variable hq. Here, 
on account of both aq and hq obeying Poisson bracket relations, Ag and fiq 
satisfy the condition. 



Xq\-q — fJ'qfJ'lq — 1- 



(57) 



However, = A_q and jjq = as will be shown below. With this change, 
using equation (57), equation (50) becomes 



bq + i 



XqiXqAq - ^.^Bq) + MA* - XqB*) 



bq + i 



Xq{XqBq- HqAq) (58) 



+Hq{lXqB* - XqA*) 



b*_q = 0. 



In order for bq to represent normal modes, the last term in equation (58) must 
vanish. This condition along with equation (57) determines A^ and /iq. After 
some straight forward algebra we find 



A, = cosh|; /xg = sinh|e^'^^ 



(59) 
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Since Aq = A_q and Bq = B_q from equations (55) and (56), it follows that 
Xq — X-q and /iq = ix-q. Substituting for Xq and ^q from equations (59) and 



(60) in equation (58) gives, 

hq + i{ujq - ir]q)bq = 0, (61) 

^l = lu, + A;r-4\Bq\']; (62) 

V,-'^iAq-A;), (63) 

with solutions 

bq (X e-^'"«*-'?«*. (64) 



In equation (64), ujq is the eigen frequency for the spin- wave normal mode and 
rjq is its damping constant. 

Next wc analyze the spin wave dynamics in the presence of oscillating terms 
in equation (50). We notice from equations (51) and (52) that terms with 
exp[±ia;t] and exp[±2ia;t] appear in equation (50). First, we suppress all other 
oscillating terms retaining only the term with exp[—2iu!t\. The equation (61) 



is modified to 

bq + l{ujq - iriq)bq + ic-'"'^' pb% = 0, (65) 

where 

P = XqiXqBl" - PqAq") + Pq{PqB'^" - XqAq"*), (66) 

i/^-^sin^^,^„^ (67) 
1 — la 
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If we assume solutions of the form bq — b°{t)e equation (65) can be 



written as 



bl + i{u;q-ujW + pb% = 0, 



(70) 



or, taking a second derivative in time, and using its complex conjugate. 



9t + (^g - ^)^ - IpI 



6^ = 0. 



(71) 



The necessary condition to see a growth of spin waves against P-mode is then 
given by 



{cug-cuf - |p|2 < 0. 



(72) 



When this condition is satisfied, the normal mode grows as 
bq oc exp (IpI^ — (uq — u})^y^^t exp[—iijUqt — r)qt\. 



(73) 



Thus, to see an exponential catastrophic growth of bq, we should have 



(74) 



With a similar argument for the term e repeating equations (65) - (74) 
we obtain a condition analogous to equation (74), 



a;, 



(75) 



with 



P = XqiXqBl' - PqAq') + Pq(PqBf - XqAq'*), 



— 27rr 

An = — 2 sin COS 6*^(1 1 ^ - l)'0o. 



1 — ia 

-27rr 

1 — ia 

27rr 



4 sin ^^cos Oqijjl, 



= — — 4sine,cosV:. 



(76) 
(77) 

(78) 

(79) 



In the above analysis, the angle 9q - the angle between the anisotropy axis and 
the x — y plane plays a crucial role. Prom equations (51) and (52), we note that 
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when 9q = 7i/2, terms with exp[±iut] vanish. Consequently, the resonance at 
uj/2 is absent. When 6q = 0, the anisotropy axis is in the x — y plane, and both 
the resonances {ug = uj/2,uj) vanish. Equations (74) with equations (66)-(69), 
and equation (75) with equations (76)-(79) provide new criteria for spin-wave 
instabilities of P-modes. The global stabihty diagrams depicted on the basis of 
this new criteria are given in Figure 5, which nicely characterize the principal 
part of the boundary curves between stable and unstable regions in Figure 3 
and 4. In particular, consistent with Figures 3 and 4 obtained from the period 
map, there is a singular structure near the vicinity of a; = (Figure 5(b'), 
5(c')) or of ha\\ — (Figure 5(b") and 5(c")) sensitive to the angle 6q. 



5 Conclusion 

We have investigated non-equilibrium states of nanoscale ferromagnets with 
uniaxial anisotropy, in the presence of an oscillating field transverse to the axis 
of anisotropy. The saturation magnetization can be driven even to the equa- 
torial plane perpendicular to the anisotropy axis. The P-modes correspond to 
new non-equilibrium states lying far from the anisotropy axis. The stability of 
the P-modes under uniform and spin wave perturbations has been studied in 
general. 

Specifically we have identified new equatorial and other fixed points, which ex- 
hibit a more complex dynamical phenomena. These novel states (fixed points) 
are realized in the case that the pumping frequency coincides with the ampli- 
tude of the static parallel field. We concentrated on the spin-wave instability 
of the above states and found that these states are always stable for the per- 
pendicular field with its strength ha± less than the value determined by the 
pumping frequency, but become unstable when ha± is "decreased" below the 
generahzed Suhl threshold. By using the period map the global stability di- 
agram in ha± — ha\\ plane is obtained, which shows an interesting singular 
structure which is very sensitive to the angle 6q between the spin-wave prop- 
agation vector (the film plane) and the direction of the uniaxial anisotropy. 

Closely following the original paper by Suhl based on Bogoliubov-type trans- 
formation, we have also obtained the generalized Suhl instability condition in 
a form of coupled spin-wave equations, which recovers the major part of the 
global phase diagram obtained by the period map. It should be noticed that 
resonance occurs at the frequencies cUg = uj/2 and cUg = cu for 9g ^ 0, but is 
absent for 9g — 0, i.e., when the axis of uniaxial anisotropy lies in the plane of 
the ferromagnet. The present findings differ qualitatively from those for the 
conventional Suhl instability near the ground state of uniaxial ferromagnets. 
Various nonlinear dynamics aspects including spin wave turbulence (see for 
example ref. [22,23]) beyond the generahzed Suhl instabihty will constitute 
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subjects which we intend to study in future. 
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Fig. 1. Stereographic projection from the unit sphere onto a complex plane. Any 
point m on the surface of the sphere can be mapped onto a point tp on the complex 
plane as shown above. Points on the upper hemisphere are mapped onto points 
inside the unit circle, while points on the lower hemisphere are mapped to the 
points outside the unit circle. All points at oo in the complex plane are identified 
with the south pole S. 

Fig. 2. Phase portrait obtained using equation (22), assuming spatial homogeneity 
of ^, for the choice h^w = uj, Keff = 0.1, D = 1, a = 0.01 and u = 0.8. The initial 

value of limjj is fixed at 0.5, while Ke^jJ is varied from —3 to 5. (a) = l + fi^, with 
fixed points corresponding to IV'oP = ^iPi = (—0.99,0.12) (unstable fixed point), 
and P2 = (0.99, -0.12) (stable fixed point), (b) 6^^ = (1 + ^^^)/2, with fixed points 
at (-.99, .12) (unstable fixed point), (0.41, -0.03) (stable fixed point), (0.99, -0.12) 
(saddle), (2.4, —0.2) (stable fixed point), (c) Enlarged image of the central portion 
in (b) showing the two fixed points, one stable, and another a saddle. In (b) and 
(c), the fixed points Pi, etc., may be identified near the respective numerical values. 

Fig. 3. Unstable regions in the (cos 6, uj) space shown by dark regions, as obtained 
from the period map M. Here, a = 0.01, Kg// = 0.1, q = 0.1. (a) sin^g = 0, 
(6) sinl?^ = 0.6, (c) sin 6*5 = 1. 

Fig. 4. Unstable regions corresponding to Figure 3 in the (/Ja||,/Ja±) space shown 
by dark regions, as obtained from the period map M. Here, a = 0.01, Kg// = 0.1, 
q = 0.1. (a) sm9q = 0, {b)sm9q = 0.6, (c) sin 0^ = 1. Not all points in the plane 
are accessible if one requires a homogeneous background, due to the condition 

Fig. 5. Unstable regions in the {cos9,u}) plane (left) and the corresponding regions 
in the {ha\i,ha±) plane (right). The parameter values correspond to the values for 
(b) and (c) in Figures 3 and 4. Instabilities at u;^ = u; are indicated by dots (•), and 
those at LOq = Lo/2 hy plus (+). For sinOq = 0, both p and p vanish as seen from 
equations (66) — (69) and (76) — (79), and hence no instability occurs. For sinOq = 1, 
the instability at Uq = uj/2 does not occur since p vanishes as seen from equations 
(76) - (79). 
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